The problem is to find the shortest distance from Bézier curve defined by the control 


points {X,Vo}.{x1.¥1},{%,¥2},{x3,¥3} to point {x,,y,}. The Bézier curve is defined in terms of 
U. 


X=Q,U+0,U+a,U+ ay 
y=b,uw+b,wW+b,utb,y 
Where a, ,d, ,4, ,dy ,b; ,b, ,b, ,and by are defined in terms of the control points. 
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The distance is s. 


s= Vix - X4) + ( - ys) 
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When the distance is minimum, the rate of change of distance is zero. 
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The order of this equation five. The solution requires trial and error methods. My 


favorite method is interpolation. Here linear interpolation will work very well if we start 
close to the correct answer. 
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Assume z is a linear function of u. 
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Now z= 0. 


The following BASIC program uses the method of finding distance. The 
program also searches for the minimum squared distance between points and 
a curve. 


REM BEZIER.BAS JIM 20DEC92 12:37 
DATA 2,3,5,8,8,14,11,17,14,17,16,15,18,11,-1 
DATA 2,10,5,12,8,11,11,8,14,6,17,5,19,10,-1 
DATA 2,5,5,7,8,8,12,12,13,14,12,17,10,18,8,17,7,14,8,12,12,8,15,7,18,5,-1 
OPEN "BEZIER.OUT" FOR OUTPUT AS #1 
OPEN "BEZ.ps" FOR OUTPUT AS #2 
CLS 
psscale = 20 
FOR example% = 1 TO 3 
REDIM rawdata(32) 
FOR I% = 0 TO 32 
READ rawdata(|%) 
IF rawdata(l%) < 0! THEN EXIT FOR 
NEXT 1% 
N% = 1% - 1 
PRINT "Example "; example%; (n% + 1) \ 2;" points" 
PRINT #1, "" 
PRINT #1, "Example "; example%; (n% + 1) \ 2; " points" 
PRINT #1,"# x y" 
J% =0 
FOR I% = 0 TO n% STEP 2 
J% =J% +1 
PRINT #1, USING "dit HH soe HHH HHH JOG: rawdata(I%); rawdata(l% + 1) 
LPRINT USING "#### #4 HHH #H# 3 0 360 arc fill"; rawdata(l%) * psscale; rawdata(I% + 1) * psscale 
PRINT #2, USING "##Ht HHH HHH HHH 3 0 360 arc fill"; rawdata(I%) * psscale; rawdata(l% + 1) * psscale 
NEXT 1% 
x0 = rawdata 
yO = rawdata 
x1 = rawdata 
y1 = rawdata' 
x2 = rawdata 
y2 = rawdata' 
x3 = rawdata(n% - 1) 
y3 = rawdata(n%) 
IF example% = 3 THEN 
‘special guess for loop 
x1=8*x1-7*x0 
yi=8*y1-7*y0 
x2 = 8* xX2-7* x3 
y2=8*y2-7*y3 
ELSE 
x1 =2*x1-x0 
yi=2*y1-y0 
X2=2*x2-Xx3 
y2=2* y2-y3 
END IF 
GOSUB distance 
LPRINT ".1 setlinewidth" 
PRINT #2, ".1 setlinewidth" 
GOSUB curveto 
e1 = totalerror 
FOR Retry% = 1 TO6 
PRINT 
PRINT "Retry "; Retry% 
PRINT #1, "Retry "; Retry% 
PRINT #1," x1 yl X2 y2 error" 
63 =.5 
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xla=x1 
DO 
x1 = x1 + (x1 - x0) * e3 
GOSUB distance 
e2 = totalerror 
IF e2 = e1 THEN 
EXIT DO 
ELSEIF e2 > e1 THEN 
x1 =xla 
63 = -e3 / 3 
IF ABS(e3) < .001 THEN EXIT DO 
ELSE 
e1 = e2 
xla=x1 
END IF 
LOOP 
e3 = 5 
yla=yl 
DO 
yi =y1 + (y1 - yO) *e3 
GOSUB distance 
e2 = totalerror 
IF e2 = e1 THEN 
EXIT DO 
ELSEIF e2 > e1 THEN 
yl=yla 
63 = -e3/ 3 
IF ABS(e3) < .01 THEN EXIT DO 
ELSE 
e1 = e2 
yla=yl 
END IF 
LOOP 
e3 = 5 
X2a = x2 
DO 
X2 = x2 + (x2 - x3) * 3 
GOSUB distance 
e2 = totalerror 
IF e2 = e1 THEN 
EXIT DO 
ELSEIF e2 > e1 THEN 
X2 = x2a 
63 = -e3/ 3 
IF ABS(e3) < .01 THEN EXIT DO 
ELSE 
e1 = e2 
X2a = x2 
END IF 
LOOP 
e3=.5 
y2a = y2 
DO 
y2 = y2 + (y2 - y3) *e3 
GOSUB distance 
e2 = totalerror 
IF e2 = e1 THEN 
EXIT DO 
ELSEIF e2 > e1 THEN 
y2 = y2a 
63 = -e3 / 3 
IF ABS(e3) < .01 THEN EXIT DO 


ELSE 
e1 = e2 
y2a = y2 
END IF 
LOOP 
IF Retry% = 6 THEN 
LPRINT "1 setlinewidth" 
PRINT #2, "1 setlinewidth" 
END IF 
GOSUB curveto 
NEXT Retry% 
LPRINT "100 200 translate" 
PRINT #2, "100 200 translate" 
NEXT example% 
LPRINT "showpage" 
PRINT #2, "showpage" 
CLOSE #1 
CLOSE #2 
END 


Bezier: 
x=a0+u* (al +u* (a2 +u* a3)) 
y=b0+u* (b1 + u * (b2 + u * b3)) 
dx4 = x - x4: dy4 = y - y4 

dx =al+u*(2*a2+u*3* a3) 
dy =b1+u*(2*b2+u*3*b3) 
z = dx * dx4 + dy * dy4 

s = dx4 * dx4 + dy4 * dy4 
RETURN 


distance: 
totalerror = 0! 
a3 = (x3 - x0 + 3 * (x1 - x2))/8 
b3 = (y3 - y0 + 3* (y1 - y2)) /8 
a2 = (x3 + x0 - x1 - x2)*3/8 
b2 = (y3 + yO - y1 - y2)* 3/8 
al = (x3 - x0) /2-a3 
b1 = (y3 - yO) /2-b3 
a0 = (x3 + x0) /2- a2 
b0 = (y3 + yO) / 2 - b2 
FOR 1% = 2 TO n% - 2 STEP 2 
x4 = rawdata(|%) 
y4 = rawdata(|% + 1) 
stepsize = 2 / (nN% + 1) 
FOR u = -1! TO 1.01 STEP stepsize 
GOSUB Bezier 
IF s = 0! THEN ul =u: z1 =z: $1 =s: EXIT FOR 
IF u=-1! THEN ul =u:z1=z:sl=s 
IFs<si THEN ul =u:z1=z:sl=s 
NEXT u 
IF s1 <> 0! THEN 
u=ul + stepsize 
IF u> 1! THEN u = 1! - stepsize 
DO 
GOSUB Bezier 
IF s = 0! THEN EXIT DO 
IF z = 0! THEN EXIT DO 
u2 =u 
72 =2Z 
temp = z2 - z1 
IF temp <> 0! THEN 
u = (z2* ul -z1 * u2)/ temp 


ELSE 
u = (ul + U2) / 2! 
END IF 
IF u> 1! THEN 
u=1! 
ELSEIF u < -1! THEN 
u=-l! 
END IF 
IF ABS(u - u2) < .0001 THEN EXIT DO 
ul =u2 
z1=22 
LOOP 
END IF 
totalerror = totalerror + s 
NEXT 1% 
PRINT totalerror; 
PRINT #1, USING "dt dH HHH HHH Ho de He He a eH x1: v1; x2: y2; totalerror 
RETURN 


curveto: 

LPRINT USING "#### ### #HH# #4## moveto"; x0 * psscale; y0 * psscale 

PRINT #2, USING "#### HHH HHHHH ### moveto"; x0 * psscale; y0 * psscale 

FG = "HHH HH HH He He He He eH HHH CUrveto stroke" 

LPRINT USING F$; x1 * psscale; y1 * psscale; x2 * psscale; y2 * psscale; x3 * psscale; y3 * psscale 
PRINT #2, USING F$; x1 * psscale; y1 * psscale; x2 * psscale; y2 * psscale; x3 * psscale; y3 * psscale 
RETURN 


